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SUMMARY 



A variational theorem is presented for a body undergoing creep. 
Solutions to problems of creep bending and creep collapse of plates, 
columns, beams, and shells can be obtained by means of the direct methods 
of the calculus of variations in conjunction with the stated theorem. 

The application of the theorem Is illustrated for plates and columns 
by the solution of two sample problems . 


INTRODUCTION 


Interest in the various effects of creep in metals has intensified 
in recent years in the aeronautical field because of actual or envisaged 
operation of aircraft in an elevated- temperature environment. Much 
experimental work has been done to determine the creep stress-strain 
relations of various metals at elevated temperatures, and much data have 
been accumulated on the creep collapse of columns and plates. Analyti- 
cally, creep problems for beams, columns, and plates are more difficult 
than the corresponding elastic problems, because the creep stress-strain 
relations (which generally include the possibility of elastic and plastic 
strains) are nonlinear. One consequence of this nonlinearity is that the 
distribution of the stresses through the thickness of the beam, plate, or 
column is not nearly so simple as in the elastic case. Mathematical anal- 
yses of these structural components have often been restricted to simpli- 
fied models (H-section columns, for example) to avoid this difficulty. 

No known systematic procedures have been devised for reducing a given 
creep problem to a set of differential equations, although some efforts 
have been made to develop variational methods for use in attacking these 
problems mathematically. Hoff (refs. 1 and 2) has suggested the use of 
a minimum-complementary-energy principle based on an analog between the 
creep-strain rates and the elastic strains in a nonlinearly elastic body; 
however, elastic and plastic strains are neglected and the method does 
not apply readily to plate and column problems. Wang and Prager (ref. 3) 
have given two general extremum principles (for a body with elastic, 
plastic, thermal, and creep strains) analogous to the minimum- complementary- 
energy and minimum-potential-energy principles in elasticity. In the 
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present paper, an analog of Reissner's variational theorem in elasticity 
(refs. 4 and 5) is formulated (for a body with elastic, plastic, and 
creep strains), and special forms of it are found to be convenient for 
handling creep problems of plates, columns, beams, and shells. The use 
of the variational theorem presented herein is illustrated by applications 
to two creep problems: creep collapse of columns and plate bending. 


SYMBOLS 


E 

I 

I l> 1 2> 1 3 

J 2 


Young's modulus 

integral to be varied (see eq. (7)) 
integrals defined by equations (46) 

second invariant of stress -de viator tensor 


M r ,M t 


Mr,0> Mt,o> M t,l 
P 

Pc 

S 

Ss 

Sd 

Ti 

Ti 

U 

U 0 


radial and tangential bending moments, respectively, 
per unit length 

radial and tangential bending-moment coefficients 
defined by equations (44) 

compressive load 

buckling load for column 

surface of a body 

part of surface where stresses are prescribed 
part of surface where displacements are prescribed 
surface traction 

prescribed surface traction on S s 
displacement of neutral surface in x-direction 
displacement coefficient in x-direction (see eqs. (28)) 


displacement of neutral surface in y-direction; also 
used as volume integral 


V 
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W 

W 0 > W 1 

R 

b 

h 

l 

m 

n j 

P 

q. 

r 

8 ij 

t 

u 


u i 


Uj: 


V 

W 

x,y,z 

7xz> ?yz* ?xy 


A 

5 i<) 


deflection of neutral surface in z-direction 
deflection coefficients (see eqs. ( 28 ) or (44)) 

radius of circular plate 

width of column 

thickness of column or plate 

length of column 

empirical constant in creep law 

unit external normal to surface of undeformed body 

empirical constant in creep law 
lateral load intensity on circular plate 
radial coordinate 
stress-deviator tensor 

time 

displacement in x-direction 
displacement vector 

prescribed displacement vector on S& 

displacement in y-direction 
displacement in z-direction 
coordinates 

shear strains in xz-, yz-, and xy-directions , 
respectively 

amplitude of initial bow of column divided by coIud 
thickness 

Kronecker delta 


e 


strain, positive in tension 
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e iJ 

€ i J 

II 

€ ij 


e r» € t 


e x> e y > e z 
9 


X 


P = r/R 
cr 


ff U 


* 

0 


* 

1 


a 0 ,cr l 


CT r* 0 t 


strain tensor 

noncreep portion of strain tensor 
creep strain tensor 

radial and tangential strains, respectively 
strain in x-, y-, and z-directions, respectively 

angular coordinate 

empirical constant in creep law 

Poisson's ratio 

stress, positive in tension 
stress tensor 

stresses defined in equation (27) 
stress coefficients in equations (28) 
radial and tangential stresses, respectively 


Dots over quantities denote differentiation with respect to time, 
a single prime denotes the elastic-plastic part of the strain, and a 
double prime denotes the creep part of the strain. 


FORMULATION OF GENERAL VARIATIONAL THEORLM 


For a body subject to creep, the equilibrium equations and strain- 
displacement relations are the same as those for an elastic body} the 
only difference in the mathematical formulation of a creep problem is in 
the stress-strain relation. In the present development, the nonlinear 
strain-displacement relations and equilibrium equations are used because 
of the intended application of the variational theorem to column and 
plate -collapse problems. However, the strains are still assumed to be 
small, and the stress-strain relations are written as for infinitesimal 



MCA IN 4003 


5 


strains. The total strain rate is separated into a creep part and a part 
which includes elastic and plastic effects. The creep part of the strain 
rate is allowed to depend on time and on the stress deviator but not on 
the stress rate or the first invariant of stress. In this paper, the 
temperature of the body is assumed to be uniform and constant in time. 

In tensor notation, a simple example of the creep law assumed is 


£ij = f(j 2 ,t)sy 


( 1 ) 


where s-y is the stress deviator and J 2 is the second invariant, 

2 s ij s ij" ^or P^P 0868 the Present paper, the only restriction on 

the stress-strain relation for the elastic-plastic part of the strain is 
that the following relation hold (see ref. 6): 


! I 

& ij d ^ij = dfir ij 


( 2 ) 


The equation for the components of finite strain is 


e 





(3) 


or, in differentiated form, 




- IK 


i + d j,i + i u k, j + u k 




The equilibrium equation is 



+ 



= 0 


(4) 


( 5 ) 
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The equation giving the surface traction T i in terms of stresses and 
displacements on the boundary is 

T ± = o Jk n^8 lk + u i>k ) (6) 

where nj is the unit external normal to the surface of the undeformed 
body* 

The components of the tensor are not ti*ue strains, and neither 

are the components of the tensor cr^ true stresses. However, the spe- 
cific virtual work due to a virtual displacement 8Ui is given by 


8W = a-y &e ii5 


Therefore, the tensor may be considered to be the generalized force 

conjugate to the generalized displacement e^. In reference 7, this 
tensor is called k^j and in reference 8 it is called cr*j. 

As in many of the variational theorems of the mathematical theory of 
plasticity (see ref. 6), the variational theorem to be stated here is for 
stress rates anfl strain rates; that is, if the Btates of stress and strain 
throughout the body are known at a given instant, then application of the 
variational theorem singles out the stress rates or strain rates that 
actually occur from those rates of stress or strain considered in the 
enunciation of the theorem. Thus, the operator & in the formal process 
of taking the variations is applied only to time derivatives of quantities 
and not to the quantities themselves. Because of this, for example, 

8e£j = 0 because depends upon the state of stress but not upon 

the stress rate. (See eq. (l).) 

The quantity to be varied was constructed by an Inverse process 
guided by the form of Reissner's theorem and the first theorem of Wang 
and Prager. By trial and error, the terms necessary to obtain the . 
desired result were found for the integral to be varied. .Let = T-j_ 
be prescribed on the part S s of the bounding surface S of a body and 
let Uj_ - u ± be prescribed on the remaining part Si of the boundary . 

The variational theorem to be proved is 81 = 0 where 



NACA TN 4003 


7 


dV - 

( 7 ) 


In "the first term, of the volume integral, e-j_ j is understood to be "written 
in terms of displacements and velocities (eq. ( 4 )). In the third term of 
the volume integral, is understood to be written in terras of stresses 

IT 

and stress rates, whereas €-j_j is written in terms of stresses. 

The variation of I is given by 



61 



8 ®iJ + ^i^k,! 6u k,,} “ 


1 

2 






dV - 



T i 8 ^i 


dS 




dS 


( 8 ) 


IT 

where 5e-y = 0 as previously noted. The second term in the volume 

integral may be transformed by an integration by parts. The following 
formula applies: 




a. .n. 


( 6 ik + 


dS - 



+ u k,i) & ij 


>1 


6u k dV 


( 9 ) 
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The third term may also be integrated by parts j thus. 


J CT i^k,i Su k,,3 dV = J ^ijpj^k, i 5 % 38 “ J ( ff i j^k, i) ^ 5u k dV ( 10 ) 


By using equations (9) , (10) , and (2), equation (8) becomes 
51 ■ 'hi ba ±i “ ^ 5 ik + u k,i) d ij] j &i k - ( a ij^k,i) ^ S *k - 

(hi + 4 ij) 5d ijj dV + J s d ij n j( 5 ik + %i) + ff ij n j\i 5 *k 


f T ± 5ut dS - f (u ± - u^STi + T ± Suj. < 

J S s ^S d L> ' -I 


or 



The coefficient of each of the variations vanishes by virtue of the 
stress-strain relations, equilibrium equations, stress boundary condi- 
tions, and displacement boundary conditions, respectively. Thus, it is 
proved that 51 = 0 for the stress rates and strain rates that actually 
occur. 


A term may be added to the volume integral to include the effects 
of body forces if necessary. In some special cases, the boundary integrals 
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should, be modified if boundary conditions other than those on Tj_ or Uj_ 
are imposed. 

In the next section, a form for I appropriate for application to 
thin-plate and column problems is developed. 


THEORETICAL APPLICATION OF VARIATIOML THEOREM 
TO ELATE AND COLUMN PROBLEMS 


The variational theorem proved in the preceding section provides a 
powerful means for deriving an appropriate two-dimensional theory for 
thin plates subject to creep (or one- dimensional theory for beams and 
columns) analogous to the two-dimensional theory for thin elastic plates. 
If certain simplifying assumptions are made about the dependence of the 
various displacements and stresses on the z- coordinate (measured normal 
to the middle surface) , the integrations with respect to z in 
the integral I may then be carried out. The equations resulting from 
51 = 0 then involve only x, y, and t as independent variables. 

In this paper, the same simplifying assumptions with regard to dis- 
placements are made that were made in the elastic case. In particular, 
the Bernoulli-Euler hypothesis is retained. This hypothesis can be 
regarded as a consequence of the equilibrium equations and the strain- 
displacement equations an applied to a thin plate. The transverse Bhear 
stresses vanish on the surface of the plate and are thus expected to be 
negligible (compared with inplane stresses) throughout the thickness. By 
any of the common stress-strain relations (elastic, plastic, or creep) , 
it follows that the transverse shear strains are negligible. If w is 
assumed to be approximately independent of z, then from 


?xz 


Bw t du 
Bx Bz 


0 


(13) 


f - 5w 5v 
y z By Bz 


0 


(14) 


it follows that 


u = U - 


Bw 

Bx 


(15) 
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v = V - 



(16) 


w = W 


(17) 


where U, V, and ¥ can now be identified as the displacements of the 
middle surface. Equations (15) , (l6), and ( 17 ) express the Bernoulli- 
Euler hypothesis. 


A more rigorous treatment, in which the dependence of the various 
quantities on the z- coordinate is deduced by expanding them in powers of 
a thickness parameter, might be possible. In some derivations of plate 
equations, U, V, and W are interpreted as weighted averages of u, 
v, and w. (See refs. 4 and 9, for example.) However, the consequences 
of assuming equations ( 15 ), (l6), and ( 17 ) to be true are expected to 
lead to accurate results for thin plates. 


Equation ( 3 ) for the finite strains in terms of displacements is 
simplified by dropping all nonlinear terms except those containing powers 
or products of slopes. This simplification leads to the Von Yi&rm&n equa- 
tions for an elastic plate. The resulting equations for the strains are 


2e-> 



dU _ 

. z 

+ l/dw') 


8x 

dx 2 

2 \8x / 


8V _ 


+ \ 

e y 

dy 

By 2 

2\by) 


_ 8u 

dV 


7xy 

dy 

H 

8x 

2z 

8x8y 


( 18 ) 


(19) 


( 20 ) 


e z ?xz “ ?yz 


0 


( 21 ) 
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In order to be consistent with these approximations, only terms conta ining 
powers or products of rates of slope are retained in the second term of 
the volume integral in I. (See eq. (T)-) 

Before the integration with respect to z can be carried out in 
the expression for I, the dependence of the stresses upon z must be 
assumed. Unlike the elastic case, the stresses do not vary linearly 
through the thickness, even though the strains do vary linearly. The 
particular advantage of the generalized form of Reissner's variational 
theorem for use in the present problem is that approximations for the 
stresses and strains can be made independently. Thus, there is no neces- 
sity to invert the stress-strain relations in order to determine appro- 
priate approximations for the stresses; in fact, to invert the stress- 
strain relations would generally be impossible. 

Many possibilities exist for assuming the form of the stresses as 
functions of z, and each possibility leads to a different set of plate 
equations; thus, no particular system of equations can be set up which 
can be called the creep equations for a plate, column, or beam. Generally, 
it is advantageous to assume the form of the stresses and displacements 
as functions of x and y also, in which case all the integrations 
in I may be performed. Setting the variation of the resulting expres- 
sion equal to zero leads to a set of ordinary differential equations in 
which the independent variable is the time. Further details of the appli- 
cation of the variational theorem to creep-col l apse and creep-bending prob- 
lems of columns and plates can best be communicated by means of examples, 
two of which are given in the following sections. 


Example of Application of Variational Theorem 
to Creep Collapse of Columns 

A simplified treatment of the problem of creep collapse of the 
uniform, pin-ended, rectangular-section column shown in figure 1 is given 
as the first example. The plastic part of the strain is omitted, and the 
stress-strain relation is assumed to be 


i = |+ Apt 5 " 1 ^ 
E 


(22) 


where the second term on the right includes the effects of primary and 
secondary creep. The axial displacement u is assumed to be given by 
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and the axial strain e by 


g _ dU _ z dfw + l/^w\ 2 
5x dx 2 


(24) 


In the present case, the following terms enter into the expression 
for I: 


• • _ •fdb , dw aV 

i ' 3€1,3 \dx 3x dx Z d x 2 


2 a iAA J = f ^ 


1 -fa w\ 2 
dx J 


> 


l(®ij + 2l ij) & ij = If + ^P tP " 1&om _ 


^ dS = -PU(l,t) = 0 (P is constant) 


(25) 


where U(0,t) = 0. The special form of I for the column Is then 


I = 


’V 


+ aw dW _ z b 2 W\ 

dx dx dx 


6x 


, + i a m 

2 2 \Bx 


_ £_ _ AptP-^o 111 
2E 


dV 


(26) 


For simplicity, the variation of a in the z-direction is assumed 
to be linear; thus, 


* z * 

a = CT n + — CT-, 
0 h 1 


( 27 ) 


The variation of displacements and stresses in the x-direction is 
assumed as follows: 
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W = hW n sin — 
U l 

U = UqX 


* _ 
0 - 



* P _ jrx 

CTn = CT-i sin — 

1 bh L l 


(28) 


where h is the thickness of the column, b is the width, nnrl l is 
the length. The dimensionless quantities on the right-hand sides of 
equations (28) are functions of time alone. The integrations in I may 
be carried out, and the result is (for m = 3) 


I = 


P 2 1 

bhE 


* I J l*0 + 3C 0 4 0 £ 


Ife 2 - i- J ■ t 


' If 0 + 5 t 01 l + 


(“o 5i o + I WA + I VVl + 6^ “l 5 ^ 


0 U 0 


(29) 


where P c is the buckling load of the column. The system of differential 

equations to which the present approximate analysis leads is obtained by 
differentiation of I with respect to each of the dotted quantities; 
thus, the result is 


ff 0 = ° 


6 *0 W 0 + % 


, P e . 

6 f w o w o 


1 * 

2 p W 0 " 


cT-j + 6a 0 W 0 — 0 

- ®o + “ & o - 7 ' ptB " lE (s) (°o 5 + I a 0°l 2 ) 

ir a i • AptP ~ lE (lr) (l “o 2 ®! + = 




0 
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The solution of this set of equations depends on the initial condi- 
tions assigned. Suppose that the column has an initial bow in the form 
of a half sine wave before the load P is applied. If the load is assumed 
to be applied rapidly, but not rapidly enough to introduce inertia effects, 
then the initial conditions for the creep problem are those existing imme- 
diately after the load is applied. If the initial bow is given by 


h A sin — 
l 


the initial conditions for the dimensionless quantities are 


w 0 (0) = — 

1 - — 

Pc 

tf 0 (0) = -1 l 


«i(°) 


12 A 



(51) 


The problem can now be reduced to solving the following equation for Wq: 




5Apt p “ 1 n (—} (w 0 + %L Wq 5 ) 

\bh ) \ ^ 5 / 


0 


(52) 


Except for differences in notation, this equation is similar to one 
given by Kempner (eq. (27) of ref. 10) which was solved by a collocation 
method for an idealized two-element column, except that the constant 27/5 
is replaced by 4/5 and p is replaced by 1. Equation (52) is easily 
integrated to give 


t p = 



w 0 2 1 + y w 0 2 (o) 

1 + ^ w 0 2 w 0 2 (0) 

5 


6 xe[— 

ibh 


2 


log 


(55) 
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The collapse time te is found by taking the limit as Wq «; thus. 


^ 1/P 



The foregoing analysis is intended merely as an example; more elab- 
orate assumptions as to the stress variation through the thick n ess and 
the variation of the several -unknowns along the length of the column lead 
to systems of nonlinear differential equations which may be solved by some 
numerical method. 


Example of Application of Variational Theorem 
to Creep Bending of Plates 

For a second example, the bending of a simply supported circular 
plate under a constant uniform load is considered. (See fig. 2.) For 
simplicity, the analysis is restricted to the time during which the deflec- 
tions are small so that the linear strain-displacement relations can apply. 
The general form for I in this case is 



Here, again, the effects of plasticity are neglected. The creep law is 
that given by equation (l). The various terms in the integrand must be 
specialized to apply to the present plate problem. In polar coordinates 
the shear stress and shear strain are zero from symmetry, and the remaining 
stresses and strains give 


^ij^ij = ^r fr r + H & t 


(36) 


If a linear strain distribution through the thickness is assumed. 


l iJ 6f iJ " 



z dW 
r dr 


& t 


(37) 
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For the elastic part of the second term. 


1 

2 


£ ±i a ±i ~ 


2E 



+ 0 + 



(38) 


For the creep part of the second term, 


hihi = f ( j 2)sijcr li3 = f(J2)j2 (39) 


Further, it is assumed that the form of f is such that 


eijffij = ^2^2 (40) 


In this example, the exponent m is taken to be unity; this assumption 
corresponds to the cubic law in the uniaxial stress case. Then, 








(41) 


The variation of the stresses through the thickness is assumed to 
be linear; thus, 


12M r z 





latfcz 


(42) 


Doubtless, after creep is well established, the stress distribution can- 
not actu ally be linear in the thickness direction, but, for purposes of 
this example, such an approximation is sufficient. When integration in the 
z-direction is carried out, 
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r 

= -f f < pM r W" + M t W* + ^Mk T 2 + M t 2 - 2^^ + 

J 0 J 0 Eh 3V / 


Z^P^M r 2 + M t 2 - M r Mt) [M r (2M r - M t ) + Mfc(2M t - M r )] (, dp d0 (4-3) 


where R is the radius of the plate, p = r/R, and the primes denote 
differentiation with respect to p. The following expressions, based on 
the elastic solution for this problem, are assumed for the displacements 
and moments: 


W = h(l - P 2)(W 0 + p 2 ^) 


M r = Eh% r?0 (l - P 2 ) { 

M t = Eh 2 ^Mfc, 0 + p3*t,l) 


(44) 


Equations (44) are substituted into equation (43) and integrated. Then, 
on differentiating with respect to the dotted quantities, the following 
five equations are obtained: 


M r ,o + 2^b,0 + Mt,l = 0 
3M r ,o + Mfc,l = 0 

f\ 


= 2 


5 It 


Wo + Wi - 2^1 jf 2M r , 0 - n(3Mt,0 + Mb,l) 

«0 - + *t,l - A.o) - 


3W 0 + Wi - - «4,o) " 6 


} (45) 
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where 


(Eh 2 ) 5 Ii = f p(l - P 2 )(2Mr - M-fc) (M r 2 - MrMt + Mt 2 )<2p 

J 0 

(Eh 2 ) 5 I 2 = f 1 ~ M r )(M r 2 - M jXt + *4t 2 )ap !> ^ 

Jo 

(Eh 2 ) 5 l 3 = f 1 P 5(a4t - Mr) (Mr 2 - Mx*4t + Mt 2 )dp 

^0 

Equations (45) can be solved for the five dotted quantities to give 


144XE 5 /R 
Wo - — 5 -(e 


Wi 


144Xe3 ( 


I 2 _ 5 C 1 + . ^ 51 3 - Ii - x 2 


|(ll - 12 + Ij) + 3(1 » ' * 2 ) 


> (4T) 


Ar,0 = Mt,0 ^(»J ' ^ - *) 


where the expressions for Iq, l2> an ^- ?5 8X6 > 

Iq = ^3M t , 0 5 - 3Mfc, 0 2M t,l “ 3M fc ,o M t,l 2 " ^A 5 ) 

I 2 = ^(W,o 5 + 9M t ,o 2 Mt,i + 7Mt,oWt,i 2 + ^A^ 

15 = i(l4M t ,o 5 + 32Mt ;0 2 Mt,q + Zfl^cPt,! 2 + ^A 5 ) 


(48) 



NACA TN 4003 


19 


Equations (47) are a system of first-order nonlinear differential equa- 
tions with time as the independent variable. The initial conditions 
required to determine the solution are given by the elastic solution to 
the loaded plate. (See ref. 11.) Thus, 


v 0 (o) = Ms * rid - a> 

0 ' E\h ] 16 


Wi(0) = - 


q/R\ 4 3(1 - |i 2 ) 
I\hy 16 


M r ,o(0) = M t ,o(0) 


> 

q./R\ 2 3 + ti 

EUi) 1 <T“ 


M t ,l(0) 


q/R \ 2 1 + 3H 

Eyh J ~l6 


J 


m 


Since M r ,0 = Mt,0 and. M r ,o(0) = Mt,o(0)> It follows that M r ^o = Mt,0* 

Equations (47) with initial conditions of equations (49) can be integrated 
numerically by the modified Euler method (ref. 12). Numerical calculations 
were made for a plate having a ratio of radius to thickness of 50. The 

value of E used was 7-4 x 10^ psi and the value of AE^ used was 
2-72 x 10° per hour. These values together with the cubic uniaxial creep 
law correspond approximately to characteristics of 2024-T3 aluminum alloy 
at a temperature of 600° F. (See ref. 2.) Results showing the time 
history of the maximum deflection of a circular plate under lateral load 
are presented in figure 3- 


CONCLUDING REMARKS 


A variational theorem for creep has been formulated which is an 
extension of a variational theorem developed by Reissner. Various systems 
of equations leading to approximate solutions to problems of the creep 
behavior of plates, columns, beams, and shells may be obtained by using 
direct methods of the calculus of variations in conjunction with the 
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stated theorem. The application of the theorem is illustrated for 
plates and columns by the solution of tvo sample problems. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 5j 1957* 
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Figure 5*- Time history of maximum deflection of simply supported circular 

TD 

plate under lateral load. — = 50 ; material, 202*t— T3 aluminum alloy at 
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6oo° F. 
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